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Abstract —In this paper, we theoretically investigate a new tech¬ 
nique for simultaneous information and power transfer (SWIPT) 
in multiple-input multiple-output (MIMO) point-to-point with 
radio frequency energy harvesting capabilities. The proposed 
technique exploits the spatial decomposition of the MIMO chan¬ 
nel and uses the eigenchannels either to convey information or 
to transfer energy. In order to generalize our study, we consider 
channel estimation error in the decomposition process and the 
interference between the eigenchannels. An optimization problem 
that minimizes the total transmitted power subject to maximum 
power per eigenchannel, information and energy constraints is 
formulated as a mixed-integer nonlinear program and solved to 
optimality using mixed-integer second-order cone programming. 
A near-optimal mixed-integer linear programming solution is also 
developed with robust computational performance. A polynomial 
complexity algorithm is further proposed for the optimal solution 
of the problem when no maximum power per eigenchannel 
constraints are imposed. In addition, a low polynomial complexity 
algorithm is developed for the power allocation problem with 
a given eigenchannel assignment, as well as a low-complexity 
heuristic for solving the eigenchannel assignment problem. 

Index Terms —RF energy harvesting, SWIPT, MIMO channel, 
singular value decomposition, channel estimation error, optimiza¬ 
tion, MISOCP, MILP. 

1. Introduction 

T he integration of renewable energy sources into com¬ 
munication networks is a hot research topic. It provides 
significant energy gains and is an efficient green communi¬ 
cation solution for the expected future data traffic increase. 
Traditional renewable energy sources such as solar energy and 
wind depend on the weather conditions and are characterized 
by high instability; the integration of these energy sources 
requires a fundamental re-design of communications systems 
in all levels of protocol stack in order to ensure robustness and 
reliability ID- Several works consider conventional renewable 
energy resources and study optimal resource allocation tech¬ 
niques for different objective functions and network topologies 
e.g., na-iQ. However, the intermittent and unpredictable 
nature of these energy sources makes energy harvesting (EH) 

S. Timotheou and I. Krikidis are with KIOS Research Center for In¬ 
telligent Systems and Networks, University of Cyprus, Cyprus (E-mail: 
{timotheou.stelios, krikidis} @ucy.ac.cy). 

S. Karachontzitis and K. Berberidis are with the Computer Engineering and 
Informatics Department, University of Patras, Greece (E-mails: {karaxodg, 
berberid} @ceid.upatras.gr) 

Part of this work was supported by the Research Promotion Foundation, 
Cyprus under the project KOYLTOYRA/BP-NE/0613/04 “Full-Duplex Radio: 
Modeling, Analysis and Design (FD-RD)”. 

This work was supported in part by the EU and national funds via the 
National Strategic Reference Framework (NSRF)-Research Funding Program 
“Thales” (project ENDECON) and in part by the University of Patras, Greece. 
Preliminary results of this work have been presented in fll . 


critical for applications where quality-of-service (QoS) is of 
paramount importance. On the other hand, recently there is a 
lot of interest to use electromagnetic radiation as a potential 
renewable energy resource. The key idea of this concept is 
that electromagnetic (EM) waves convey energy, which can 
be converted to DC-voltage by using specific rectenna circuits 

ii, B- 

Most of the work on radio-frequency (RE) energy transfer 
concerns the design of rectenna circuits in different frequency 
bands, which is a fundamental aspect towards the development 
of this technology ifTOt - lfT^ . Erom an information theoretic 
standpoint, evaluating the channel capacity for different net¬ 
work configurations, when RE EH requirements characterize 
the receiver nodes, is a challenging problem. The work in 
|[T^ discusses the joint transfer of information and energy for 
a single-input single-output (SISO) channel and is extended in 
for a set of parallel point-to-point channels; the authors 
in ca study the capacity for two baseline multi-user systems 
with RE energy constraints, namely multiple access and mul¬ 
tihop channels. However, information theoretic studies assume 
that the receivers are able to decode information and harvest 
energy from the same RE signal without limitations. Although 
this assumption provides some useful theoretical bounds, it 
cannot be supported by the current practical implementations. 

In order to satisfy the above practical limitation, the work in 
|[T^ deals with the beamforming design for a basic broadcast 
multiple-input multiple-output (MIMO) channel, where the 
source conveys information to one receiver and transfers 
energy to the other one. In that work, the authors introduce 
two main practical techniques for simultaneous information 
and power transfer (SWIPT): a) “time switching” (TS), where 
the receiver switches between decoding information and har¬ 
vesting energy and b) “power splitting” (PS), where the 
receiver splits the received RE signal in two parts for decoding 
information and harvesting energy, respectively. This work 
is extended in C3 for imperfect channel information at the 
transmitter using robust optimization. 

The employment of the above two practical approaches in 
different fundamental network structures, is a hot research 
topic and several recent works appear in the literature. The 
works in di-Eil focus on the TS technique for different 
network topologies. Specifically, in |[T^ the authors investigate 
the optimal switching strategy for a SISO channel in order 
to achieve various trade-offs between wireless information 
transfer and EH with/without channel state information (CSI) 
knowledge at the transmitter. The work in CD presents a time 
division multiple access (TDMA)-based multiuser broadcast 
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channel where the downlink (broadcast channel) is used for 
EH while the uplink for conveying information from the users 
to the access point. In 1^ . the authors apply the TS technique 
for a basic relay channel and investigate the optimal switching 
rule. A TS switching architecture is proposed in (Ell for a 
multiple-input single-output (MISO) system, where the single¬ 
antenna receiver uses its harvested energy to feed a quantized 
CSI. On the other hand, in ll22l the authors study the perfor¬ 
mance of a cooperative system, where the relay is powered by 
employing a PS technique on the received signal. The work 
in ll^ studies the optimal transmitted power for a MISO 
interference channel with PS where the destinations have both 
information/energy constraints. The TS and PS techniques 
have also been exploited in cognitive radio networks to achieve 
joint information and energy cooperation between a primary 
and a secondary system, resulting in improved performance 
for both systems due to the added RF EH capabilities 1^ . 
Recently, the authors in ll^ take into account the energy 
consumption of the active electronic circuits (e.g., mixers) 
and they propose an integrated receiver architecture, where 
PS is employed at the baseband domain. A recent overview 
of SWIPT systems with a particular focus on the hardware 
realization of rectenna circuits and practical techniques that 
achieve SWIPT in the domains of time, power and antennas 
is provided in ll^ . 

Although most of the previous studies refer to specific and 
fix network topologies, more recent works study TS/PS for 
large scale networks with location randomness. In llZ7]| . the 
authors study the interaction between primary and cognitive 
radio networks, where cognitive radio nodes can harvest 
energy from the primary transmissions, by modeling node 
locations as Poisson point processes. A cooperative network 
with multiple transmitter-receiver pairs and a single EH relay 
is studied in 1^ by taking the spatial randomness of user 
locations into consideration. In ||29l, the authors study the 
performance of a large scale ad hoc network with/without 
relaying where receivers have PS capabilities. 

In contrast to the conventional TS and PS approaches, we 
propose a new technique for SWIPT in the spatial domain for 
a basic point-to-point MIMO channel. In this work, spatial 
domain does not refer to the antenna elements ca, oa but 
mainly on the spatial degrees of freedom of the channel. Based 
on the singular value decomposition (SVD) of the MIMO 
channel, the communication link is transformed to parallel 
channels that can convey either information data or energy; 
this binary allocation is in respect to the current practical 
limitations. In order to make our analysis more general (and 
practical), we assume an imperfect channel estimation that 
affects the orthogonality of the eigenchannels (the imperfect 
channel knowledge generates an interference to the eigenchan¬ 
nels). We study the minimization of the transmitted power 
when the receiver is characterized by maximum power per 
eigenchannel, information rate and power transfer constraints 
which is a mixed-integer nonlinear optimization problem and 
propose several solution methods. It is worth noting that the 
main purpose of this work is to introduce a new technique for 
SWIPT in MIMO systems; this technique is studied from a 
theoretical standpoint and practical implementation is beyond 


the scope of this paper. 

More specifically, the contributions of this paper are: 

• Construction of mixed-integer second-order cone pro¬ 
graming (MISOCP) and mixed-integer linear program¬ 
ming (MILP) formulations to the general problem that 
can be solved using standard optimization tools to provide 
optimal and near-optimal solutions, respectively. 

• Development of an optimal polynomial complexity algo¬ 
rithm for the solution of the special problem of having 
no maximum power per eigenchannel constraints. A low 
complexity polynomial algorithm is also developed for 
the general problem with excellent performance. 

• Derivation of a waterfilling-like procedure for the solu¬ 
tion of the optimal power allocation problem with fixed 
eigenchannel assignment, of low polynomial complexity. 

Notation: Upper case and lower case bold symbols denote 
matrices and vectors, respectively, while calligraphic letters 
denote sets, diag(x) represents a diagonal matrix with the 
vector X in the main diagonal, det(-) denotes determinant. 
In is the identity matrix of order n, log(-) denotes the 
logarithm of base 2, E['] represents the expectation operator, 
the superscripts H and * denote the Hermitian transpose and 
conjugate operations, respectively. \A\ denotes the cardinality 
of set A, while ^ \ S is the difference between sets A and B. 

The paper is organized as follows. Section II presents the 
system model, introduces the spatial domain EH and the 
associated optimization problem. Section III deals with the 
optimal solution of (a) the general problem considered using 
MISOCP and MILP formulations and (b) the special case with 
no maximum power per eigenchannel. Section IV develops 
an optimal procedure for the power allocation problem for a 
given eigenchannel assignment. Section V introduces a low- 
complexity suboptimal heuristic for the considered optimiza¬ 
tion problem. Section VI discusses the numerical performance 
of the proposed schemes and Section VII concludes the paper. 

H. System model & problem formulation 

We assume a simple point-to-point MIMO model consisting 
of one source with Nt transmit antennas and one desti¬ 
nation with Nr receive antennas. The source is connected 
to a constant power supply while the destination has RF 
transfer capabilities and can harvest energy from the received 
electromagnetic radiation. We consider a fiat fading spatially 
uncorrelated Rayleigh MIMO channel where H G 
denotes the channel matrix. The channel remains constant 
during one transmission time-slot and changes independently 
from one slot to the next. The entries of H are assumed to be 
independent, zero-mean circularly symmetric complex Gaus¬ 
sian (ZMCSCG) random variables with unit variance (which 
ensures rank(H) = N = min{A^T, Nr})- The received signal 
is described by 

y = Hx + n, (1) 

where x G denotes the transmitted signal with 

E[xx^] = Q and n G represents the noise vector 

having ZMCSCG entries of unit variance. We assume that the 
channel matrix is subject to a channel estimation error and 
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Fig. 1. The SVD transformation of the MIMO channel into N parallel 
AWGN channels. 

therefore is imperfectly known at both the transmitter and the 
receiver with a MMSE estimation error E = H — H, where 
the entries of E are ZMCSCG with variance and the entries 
of H are also independent and identically distributed (i.i.d.) 
ZMCSCG with variance l-cr^ llITIl . The channel estimation is 
performed at the destination via a downlink pilot sequence and 
is communicated to the source by using an uplink feedback 
channel; the channel estimation process is beyond the scope of 
this paper. We note that the parameter cr‘1 captures the quality 
of the channel estimation and is assumed to be known to both 
the transmitter and the receiver. The work in ED provides 
some examples for related to different channel estimation 
schemes and channel statistics; in this paper, we assume that 
cTg is a constant without further implications. Based on the 
estimated channel H, a lower bound of the instantaneous 
mutual information is given by ED 

/(x; y) = log det . (2) 

It is worth noting that for = 0 (perfect channel knowledge), 
the above expression gives the exact mutual information of 
the MIMO channel. By using the SVD of the H channel, it 
has been proven in ll3T]| that the lower bound of the mutual 
information is maximized when Q = diag(Pi,..., Pnt) 
takes the form 

/(x;y) = 5:iog(l + ^^^), (3) 

zGA/ fc / 

where \i,i E Af = {I,-.-, N}, is the i-th eigenvalue of HH^, 
Pi is the power allocated to the i-th eigenchannel and P = 
Pi- It is further assumed that Ai > A 2 > ... > Xn- 
This bound can be achieved using SVD on the estimated 
matrix H and by appropriate precoding and receiver shaping 
at the source and the destination, respectively. Fig. [T] schemat¬ 
ically presents the system model and the transformation of the 
MIMO channel to N eigenchannels for potential transfer of 
information/energy. We show the SVD coordinate transforma¬ 
tions when the matrix H has a SVD equal to H = UEV* , 
where U G and V G are unitary matrices 

and E G is a diagonal matrix with elements the 

singular values of the matrix H i.e., ai, i E Af and ai = 
Based on the SVD, the channel input x comes from the 
transformation of the original sequence x with the matrix V 
while the received sequence y is transformed to the sequence 
y by using the matrix U. In the ideal case where = 0, the 
MIMO channel is decomposed to N parallel SISO channels 


with 

Pi = aiXi -h fii^ with i E Af^ (4) 

where fii denotes the AWGN for the i-th parallel channel 
and has the same distribution with rii (due to the unitary 
transformation) in reception. 

The expression in ([3]) shows that the imperfect channel 
estimation generates an interference to the N parallel SISO 
channels (eigenchannels). We assume that the destination is 
characterized by both information rate and RF EH require¬ 
ments; this means that for the duration of one transmission 
the destination requires an information rate Cj and energy 
Ceh = Ceh 15 as input to its rectenna circuits (the amount 
of energy that can be stored, Ceh, depends on the energy 
conversion efficiency of the specific implementation, J) e.g., 
1^ . (321. For simplicity and without loss of generality we 
assume in the mathematical analysis that S = 1. 

A. SWIPT optimization problem 

The proposed scheme exploits the SVD structure of the 
MIMO channel and achieves SWIPT in the spatial domain. 
More specifically, the transformation of the MIMO channel to 
N parallel SISO channels allows the simultaneous transfer of 
data traffic and RF energy by using an eigenchannel either to 
convey information or energy. An eigenchannel cannot be used 
to convey both information and energy; this limitation refers 
to practical constraints and is inline with the other approaches 
proposed in the literature, e.g. power splitting. Based on the 
diagram in Fig. [B at the output of each eigenchannel there 
is a switcher which drives the channel output either to the 
conventional decoding circuit or to the rectification circuit. 
During each transmission an appropriate optimization problem 
is solved, which determines the usage of each antenna and a 
switching mechanism selects the appropriate circuits. 

In this paper, we focus on the minimization of the transmit¬ 
ted power given that both information and energy constraints 
are satisfied. Based on the notation considered, the proposed 
technique introduces the following optimization problem 


min P = Pi (5a) 

ieV 

s.t. + >C,, (5b) 

ieA/ ^ fc / 

^(1 - + alP) > Ceh, ( 5 c ) 

ieM 

^ ^ Pi ^ Pmax-j G{0,1}, ieM, (5d) 


where Pmax indicates the maximum power that can be used 
in each eigenchannel, while binary variable tt^ indicates 
whether the i-th eigenchannel is used for information transfer 
{Ei = 1) or energy transfer (tt^ = 0). Note that terms 
log equal to the more intuitive representation 

Ei log -f when Ei E {0,1}, i E Af. This mathe¬ 

matical program involves binary and continuous variables, as 
well as nonlinear functions; hence it belongs to the class of 
mixed-integer nonlinear optimization problems, which are very 
hard to solve in general. 
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The incorporation of into ^ makes the problem more 
realistic and general, with interesting implication effects. On 
the one hand, increasing reduces the achieved information 
rate in dSB. On the other hand, the added cross-eigenchannel 
interference becomes a useful source of energy transfer in 
(l5cl) . which implies that even zero-power EH eigenchannels 
can contribute to the satisfaction of the EH constraint. Erom an 
optimization point of view, the problem becomes more difficult 
to tackle due to the presence of term cr^P in the denominator 
of the logarithmic terms in (l5bl) . Hence, the developed solution 
methodologies are quite general and can be applied not only 
for the solution of the special case with = 0 , but also for 
other problems with similar information rate expressions. 

We note that the energy harvested due to the receiver noise 
is negligible. In addition, the optimization problem can be 
solved at the source node and then control bits can be used 
in order to inform the receiver about the content of each 
eigenchannel. 

B. Implementation issues 

The main challenge to implement the proposed SWIFT 
technique is to perform the required signal processing, in¬ 
cluding the eigenchannel decomposition, with ultra-low power 
consumption. This requires the use of analog passive electronic 
elements (e.g., Schottky diode), which is in line with current 
implementations of conventional wireless power transfer ||3^ . 
El. In respect to these requirements, analog eigenmode 
beamforming with passive electronic elements (180-degree 
hybrid couplers) is a promising technology recently proposed 
ES, that could be used for the practical implementation 
of the proposed technique. This implementation technology 
achieves diagonalization of the MIMO channel matrix without 
destroying the energy content of the received signal. 

It is worth noting that although digital beamforming dom¬ 
inates current communication systems, analog beamforming 
is an emerging technology with particular impact on future 
millimeter-wave (mmWave)- based systems Ea, El- In these 
systems, the transceivers are composed of a large number 
of antennas and the conventional baseband/digital processing 
requires a large number of analog-to-digital converters and 
RE chains, which results in a high complexity and power 
consumption. Also note that SWIFT is only studied theo¬ 
retically and not used in current devices; proof-of-concept 
implementations are still under investigation. 

HI. Optimal Eigenchannel Assignment and Fower 
Allocation 

The considered problem is nonlinear and combinatorial in 
nature, and hence very hard to solve. Moreover, standard 
MILF/MISOCF solvers are not directly applicable for its 
solution due to the presence of nonlinear terms. In this section, 
two formulations are proposed for the optimal/near-optimal 
solution of problem ^ using standard solvers. The first, 
transforms the problem into an equivalent MISOCF formu¬ 
lation that provides the optimal solution, while the second 
approximates the logarithmic functions using piecewise linear 
approximation (FLA) to arrive at a near-optimal (in most cases 
optimal) result by solving a series of MILF problems. These 


approaches have exponential complexity and are only used for 
benchmarking purposes; the reason for discussing both is that 
the MISOCF guarantees optimality, while the MILF has more 
robust performance (see Section [Vl]). In addition, we provide 
a polynomial algorithm that optimally solves the problem for 
the special case of Pmax = oo. 

A. Optimal MISOCP solution 

Towards the optimal solution of (|5]), we write it as: 


min P = Pi 

( 6 a) 

ieAf 



( 6 b) 

iEAf ^ e / 


y^iPi - nPi)K + NalP - al Y. TT, 

IV 

ieM ieAf 

( 6 c) 

0<Pi< Pmax, n e {0, 1} , ie N. 

( 6 d) 


Formulation ([ 6 l) is not MISOCF due to the products of 
variables, i^iPi and tt^P, i G A/", in ([ 6 bl)-([ 6 cb and the general 
product of terms in ([ 6 bb . Next we show how these can be 
transformed into equivalent representations of MISOCF form. 
Firstly, we state two auxiliary propositions. 

Proposition 1. Let w G {0,1} and x G 
Constraint y = wx can be equivalently represented with the 
following linear constraints: 

wx^^ <y<x — {l — w)x^^ ^ (7a) 

X — {1 — w)x^^ < ^ < wx^^, (7b) 

Proof: The proof is easily derived by checking that w = 0 
and w = 1 yield ^ = 0 and y = x, respectively. ■ 

Let us define Pf = i^iPi and P/^^ = ttiP, i E ff, which 
appear in constraints (Ibbl) and (Ibcl) . Based on proposition [B 
equalities Pf = TiiPi, i e ff can be represented by constraints 
(iTab-dTbb with y = Pf, w = iTi, x = Pi e [0, Pmax]', similarly 
we can represent P/^^ = TTiP, with y = P/^^ w = iTi, 
X = P e [0, NPmax]- Next, we will show that the general 
product appearing in ([ 6 bb can be represented by a hierarchy of 
convex second-order rotated cone (SORC) constraints of the 
form X 1 X 2 > Xq, > 0 , X 2 > 0 using proposition [ 2 l 

Proposition 2. ( 051/ . p.l05): The geometric mean constraint 
(GMC) X 1 ...X 21 , Xi >0, i = 1,..., 2^ is convex and can 

be represented by a hierarchy of SORC constraints. Let us 
define existing variables as Xi = xo^i and new variables Xk^i > 
0, /c = 1,...,/, i = Then, the GMC is equivalent 

to the following SORC constraints 

layer k: Xk-i,2i-iXk-i,2i > = 1 ,fc = 

xi,i > t. 

The proof is based on the fact that each of these constraints 
is SORC, hence convex, while if these constraints hold the 
GMC holds as well. To see why this is true, let us consider 
the case X 1 X 2 X 3 X 4 > Setting Xi = xo, 2 , i = this 
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constraint becomes: 


^0,1^0,2 >^1,1, ^0,3^0,4 > (8a) 

3^1,12^1.2 >2^2.1, 2^2,1 > t. (8b) 

The two constraints in (ISat yield that xq, 1 X 0 , 2 ^ 0 , 3 ^ 0,4 > 
1^1 2 ’ which combined with (ISbl) gives x\ix\ 2 ^^21 ^ 
hence, if (ISat-dSbl) hold, the original constraint also holds. 
To bring (l6bl) into the GMC form, we define M = [log A^] 
and Nu = 2^, and rewrite it as: 

N Nu-N 

n(l+(7?P + P,^A,) H 

i=l i=l 

Following proposition [21 we can represent (l6bl) with the hierar¬ 
chy of SORC constraints by setting xo,i = lpa‘^PpP-^Xi > 
0, i e Af, xo,i = 1 + cTgP > 0, i = + 1, Nu and t = 

2 CiINu (i+cr2p) Towards this direction, a total of A^^ —1 new 
variables are needed, Xk,i > O^k = i = 1 ,..., 2^“^, 

and Nu/2 SORC constraints. 

Based on the above analysis, the original problem can be 
rewritten as: 


min P = Pi (9a) 

ieN 

Xk-i,2i-iXk-i,2i>x‘l^i, A; = 1,..., M,z = 1,..., 2^“^, (9b) 
xm,i > t, t = + cr^P), (9c) 

xo,i = 1 + cTgP + Pj Ai, i G A/*, (9d) 

xo,z = 1 + cTgP, z = A^ + 1,..., A^^,, (9e) 

> 0, /c = 0, ..., M, z = 1, ..., 2'^-^ (9f) 

^(P. - p:)\i + (7,"(iVP - ^ Pt^) > Ceh, (9g) 

iGAA ieM 

Pi < TTiPmax, P^ < Pi, (9h) 

P: > Pi -{I- TTi)P,nax, > 0, i € V, (9i) 

Pt^ < TiiNPmax, P‘°‘ <P, ie AT, (9j) 

ptot > p _ _ 7ri)iVP^,„ p/°* > 0, i € V, (9k) 

0 ^ Pi ^ Pmax-i G {0,1} , z G A/”. (91) 


In formulation (|9]), (l9bl)-(l9fl) represent the information con¬ 
straint, ( |9g| the EH constraint, while (l9hl) . (l9il) and (|^, (|9B 
are used to linearize equations P^ = Tr^P^ and P/^^ = tt^P, 
respectively. Note that formulation (|9]) is MISOCP and can be 
solved with standard optimization solvers. 


B. Near-optimal MILP solution 

In this section we propose an iterative approach for the 
solution of (O, where in each iteration a MILP problem is 
solved. Although, this approach yields near-optimal results (in 
most cases optimal) it is considerably faster than MISOCP 
for large problems. Towards this direction, we examine the 
implications of having a fixed value for the total power P in 
(l5bl) . Assuming that this value is equal to P/, the problem 



Fig. 2. Piecewise linear approximation; = 2; and \i/{l-\- cr’^Pf) = 

10 . 


becomes: 


min P = Pi 

(lOa) 

ieN 



( 10 b) 

y^(l - TTi) {PiXi + cr‘^P) > Ceh, 

( 10 c) 

ieN 


0 ^ < Pmaxi ^ {O 5 1 }: '1 ^ Af. 

(lOd) 

Using the transformations, P[ = tt^P^ and P^^^ 

= TTjP, 


introduced in section IIII-AI via proposition 1, all the con¬ 
straints of problem ([TOb become of MILP form, except from 
constraint (fTOhli . due to the presence of terms fi{Pi) = 
log ^1 + i^a^Pf ) • Nonetheless, these terms can be trans¬ 
formed into MILP form by considering appropriate piecewise 
linear approximation (PL A) functions gi{Pi) ~ fi{P-^). PL A 
functions gi{P^) can also be selected to provide either upper 
or lower bounds (UB or LB) to the optimal value of problem 

m, if giiPi^) < MP) or gi{Pn > ppp 0<Pr< 

Pmax, 'I ^ A/", respectively. Constraints gi{Pi ) < MP) and 
gi{Pi ) ^ fi{Pi ) be achieved by taking the lower envelop 
of a number of lines that intersect fi{Pi) from below or are 
tangent to fi{Pi) from above respectively, as illustrated in 
Fig. El 

One can easily find the PLA of a logarithmic function, by 
considering a fixed number of points, at integer multiples of a 
parameter n, such that segments are either tangent (for UB) or 
intersect (for LB) the logarithmic function. Nonetheless, this 
approach has no control over the approximation error, while 
the maximum error for different segments varies. For this 
reason, an approach has been followed that produces segments 
when needed so that the approximation error does not exceed 
a predefined value emax lESI- Although this approach may 
lead to a large number of segments we can adjust emax to 
achieve the desirable number of segments and at the same 
time maintain the error less than a constant known value. 
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The approximation procedure, involves two steps for the 
production of each segment. In LB-PLA, where an LB en¬ 
velope needs to be produced, the first step for the deriva¬ 
tion of segment /, starts from a known intersection point 
{Pill J(Pi",i)} (initially, = 0}) and finds a 

line ai^iP^ +/3i^i such that< emaxJor 
f{Pi) ^ ^i,iPi + Pi,I' the second step, a new intersection 
point is found, such that = (yi,iP^^i^i + Pi,u which 

completes the creation of the Ith segment. In UB-PLA, where 
a UB envelope needs to be produced, the first step for the 
derivation of segment /, starts from a point of maximum 
error, + emax} (initially, the procedure starts 

from {P[i = 0,/(P/^i) = 0} to produce such a point) and 
computes a line OLi^iP^ -f bi^i that is tangent to /(P/^). In the 
second step, a new point of maximum error is found such that 
^i,iPP,i+i + Pi,I = fiPP,i+i) + ^rnax, which Completes the 
creation of the Ith segment. This two-step iterative procedure 
is followed until Pmax is reached. 


Having obtained PLA functions gi{P^), i G A/', we now 
describe how to obtain a problem formulation with MILP 
constraints. Let us assume that function fi{P^) is approx¬ 
imated by Li linear segments with slopes and 

start-points P/^^,...,P/^^.. Let us also assume that P^l.^i = 
Pmax- Because fi{Pi) are monotonically increasing and 
concave, the PLA functions gi{P^), will also be concave, 
with segments whose slopes have monotone decreasing values: 
(^i,i > Q^i, 2 > ••• >Oii,Li- Let = l,...,Pi be the value 
of P[ corresponding to the Ith linear segment so that 0 < 
^i,i < Pu-\-i ~ Piu ^ ~ 1, Under the assumption that 

ii,m = Plm+ 1 -Plm^^ = 1 when > 0, it is true 

that ff = and also that gi{Pi) = 

In other words, P^ can be replaced by the sum of variables 
= 1,...,P^ if we can ensure that the solution of the 
optimization problem will always be such that each is 
nonzero only when the variables m = — 1 have 

obtained their maximum value. 


To explain why this condition holds for problem cs, 
note that variables P^ should be as small as possible to 
yield minimum total power, while the impact of gi{P^), and 
hence of variables should be as large as possible to 
ensure satisfaction of constraint (llObb with minimum power. 
Because the slopes of the PLA segments have monotonically 
decreasing values, the first segment has the largest impact; 
this implies that will be the first variable to be assigned a 
nonzero value and only if becomes equal to its maximum 
value, variable ^^,2 will be assigned a nonzero value. This 
ensures the assumption made that will be non-zero only 
if = 1, •••,/ — 1 have attained their maximum value. 


Algorithm 1 : MILP solution to problem (|5]) 

1: Init. Plb = 0; Pub = NPmax\ e = 10“^; Cmax = 10“^; 

/ = o. 

2: while \PuB - Plb\ > e do 
3: Set Pf = {Pub + Plb)I‘^- 

4: Compute aij and P^i, i e JV, I = 1, ...,1/^ according 

to Cmax using PLA. 

5: Solve problem (fTTI) using a MILP solver to obtain Pi, 

TTi, i e JV. 

6: if (problem feasible) AND (P < Pf) then 

7: Set Pub = Pf \ tt* = tt^; / = 1. 

8 : else 

9: Set Pub = Pf- 

10 : if / = 1 then 

11: Solve the power allocation problem with fixed assign¬ 

ment TT^, i e JV (Section HVl). 

12: else 

13: Problem (|5]) is deemed infeasible. 


Based on the above discussion formulation (fTOl) becomes: 

min P = E Pi (11a) 

ieAf 

Li 

s.t. > Cl, (11b) 

ieAf 1=1 

Li 

PI'= iefC, (lie) 

1=1 

0 < ^i,i < PI'^i+i - PU ieN, 1 = 1 ,..., Li, (lid) 

Constraints @-®. (lie) 

Note that formulation (UB is MILP, as constraints dllal) - 
(fTTdli . which represent constraint (II Obi) , also become linear 
under the considered PLA scheme. Also, by considering UB- 
PLA (LB-PLA) for fi{P^) a lower-bound (upper bound) to 
problem ([Tob is obtained, as we need to use more (less) power 
to satisfy the information constraint. 

If the solution of formulation STT\f is feasible and the total 
required power is less than Pf, i.e. P < Pf, then Pf is 
an upper bound to the optimal total power of problem ([5]), 
otherwise Pf is a lower bound. Based on this observation, we 
can construct a bisection procedure to obtain a near-optimal 
value for the total power, according to algorithm [T] 

Loss of optimality in algorithm [T] is primarily due to the 
PLA of the logarithmic terms in the total information rate 
expression which may result in a non-optimal eigenchannel 
assignment. Note that global optimality is preserved when 
the optimal eigenchannel assignment is found, by optimally 
solving the power allocation problem with fixed eigenchannel 
assignment (line 11). 

C. Special Case: Pmax = oc 

In this section, we deal with the solution of problem ([5]) 
when Pmax = CO. We show that we can obtain the optimal 
assignment by examining a polynomial number of information 
and EH eigenchannel combinations and develop an algorithm 
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that solves this problem optimally in 0{N‘^). We also derive 
analytical results for obtaining the optimal power allocation 
for a fixed eigenchannel assignment, as well as for deciding 
whether a candidate assignment is potentially optimal prior to 
computing the power allocation. 

Assume that a given assignment of channels is made such 


that TTi^ = 1, ii e X and 77^2 = 0, 22 ^ where X and £ 
denote the sets of information and EH eigenchannels, respec¬ 
tively. Due to this assignment, problem (O, with Pmax = co, 
can be written as: 

min P = Pi (12a) 

ieAf 

s.t.^ log (1 + > Cl, (12b) 

hex ^ 

(^ 22^22 P CeH, ( 12 c ) 

he^ 

Pi > 0, i G AT. (12d) 


Three important remarks can be made regarding (fT^ . 
Remark 1. At most one EH eigenchannel can have nonzero 
power. 

Remark 1 is true because there is no imposed upper bound 
on the transmitted power in each eigenchannel. As a result, 
if a set of eigenchannels £ are assigned to satisfy the EH 
constraint, only the one with the largest eigenvalue in this 
set, Ae, e G f will be nonzero. Despite the fact that the 
power assigned to other channels in £ is zero, they still 
contribute towards the satisfaction of the EH constraint due 
to the presence of the term g\P in the EH constraint. 
Remark 2. Eigenchannels with P^ = 0 are those with the 
overall smallest eigenvalues. 

This is because “better” eigenchannels are beneficial for 
both the satisfaction of the information and EH constraints. 
Hence, if we consider two eigenchannel with > Aj, P^ = 0 
and Pj = c > 0 , by simply switching the power of the two 
eigenchannel, i.e., Pj = 0 and P^ = c > 0, more benefit can 
be obtaining for the corresponding constraint without affecting 
the other constraint or the total power in any way. 

Remark 3. At the optimal solution, no information eigenchan¬ 
nel should have zero power. 

To see why this is true, notice that eigenchannels with zero 
power can contribute to the EH constraint but not to the infor¬ 
mation constraint. This is a result of the cross eigenchannel 
interference, cr^P, that is present in all eigenchannels due to 
the channel estimation error. Hence, zero power information 
eigenchannels should always be used as EH eigenchannels to 
have a positive contribution towards the satisfaction of the EH 
constraint. 

These remarks imply that in order to find the optimal 
assignment we need two different indices. The first index e 
indicates the energy eigenchannel with a possibly non-zero 
value, while the second indicates the channel with the largest 
eigenvalue such that Pi^ = 0, 22 ^ ^ \ Based on the 
above observations, we can conclude that in order to find the 
optimal assignment, we only need to examine 0{N‘^) different 
assignment combinations. The next result, summarizes the 
solution of each assignment combination. 


Theorem 1. Let X, £ be a fixed assignment, e = min{f}, 

C = - W Ti = {i:iel, Pi> 0}, 

and /3 be given by: 


(3= i2^^/'[[Xi 


WT] 


(13) 


ieXi 


The optimal total power for this assignment is given by P = 
{y — y > 1, where: 


y = 


(G£;jj' + |£:|)/Ae+0'^ ^ 

C+l^l/^e 


*/<Trvc>^ + i 

otherwise. 


(14) 


In addition, given P and y, the optimal power allocation for 
a given assignment, X, £, is: 


Pi = max{0,^(/3 - l/\i)}fi G X, 



7 = e, 
i G £ \ {e} 


(15) 

(16) 


Proof: The proof can be derived by following similar 
analysis with Section III in (U. ■ 

Theorem 1 consists of two parts; the first part illustrates 
how the optimal power, P, of the examined assignment can 
be obtained, while the second part provides its optimal power 
allocation. To obtain P it is essential to compute parameter 
f3 which requires finding set Xi containing the indices of 
information eigenchannels with non-zero power. These can be 
obtained by sequentially adding information eigenchannels to 
/3, starting from the best one, until an eigenchannel, k, is found 
for which /3 — 1/Xk < 0, k e X. One important observation 
regarding theorem 1, which emanates from remark 3, is that 
condition P — IjXj > 0, j = max{X} is necessary for a 
candidate assignment to be optimal. Hence, any assignments 
that do not satisfy this condition are immediately rejected 
without any further consideration. 

To check for infeasible solutions, we further need to ensure 
that the obtained solution for y yields a positive value for Pe 
and that ^ > 1. As ^ is a monotonically increasing function of 
P, by finding the best feasible value for y among all feasible 
assignments, yopt, we can derive the optimal power allocation 
of ([T^ . as outlined in algorithm [2l 

It should be emphasized that algorithm [2] has two very at¬ 
tractive characteristics: (a) it solves a nonlinear combinatorial 
optimization problem involving binary variables, in polyno¬ 
mial time, as it requires the examination of a polynomial 
number of fixed assignments (approximately N‘^/2), and (b) 
the optimal power allocation needs to be derived only for 
the optimal assignment Xopt, eopt, and not for all examined 
fixed assignments which reduces the execution time of the 
algorithm. Although each assignment appears to be of com¬ 
putational complexity 0{N) due to the presence of 
and riiex Xi, we can reduce the computational complexity to 
0(1). This can be achieved by storing the sum and product 
terms for fixed EH assignment, e, and updating their values for 
an increasing number of information channels i = 6 + 1,...,^". 
Hence, the total complexity of algorithm [2] is 0{N‘^). 





Algorithm 2 : Optimal solution to problem ([T^ 

1: Init. yopt = oo, Xopt = 0, and Copt = 0. 

2: for e = 1 to do 

3: for i = e + 1 to do 

4: Initialise fixed assignment: X= — 

5: if (|X| > 0) then 

6: Compute f3 according to (fTSl) . 

7 : if (/3 — 1/Xj > 0, j = max{X}) then 

8: Compute the value of y according to (O. 

9: if ({yopt > y) and (y > 1)) then 

10: Store the optimal solution found so far: 

11: Vopt — Xopt — X, and ^opt ~ 

12: if (1 < yopt < oo) then 

13: Having found optimal assignment Xopt, ^opt and yopt 

compute optimal power allocation according to (fTSl) and 

(HSl). 

14: else 

15: Deem problem infeasible. 


Algorithm 3 : Optimal power allocation for problem ([TTb 

1: Init. Plb = 0; Pub = NPmax\ e = 10“^; / = 0. 

2: while \PuB - Plb\ > e do 
3: Set Pf = (Pub + Plb)/‘^- 

4: Solve the information subproblem (fTSl) according to the 

procedure discussed in section IIV-AI 
5: Find the optimal power allocations, Pi,i G {X U f} 

according to cases 1-4 discussed in section HV^ 

6: if (problem feasible)AND(^.^|^^^| Pi < Pf) then 

7 : StiPuB=PpPt =Pi\ f = X 

8 : else 

9 : Set Pub = Pf- 

10: if / = 0 then 

11: Problem dnii is infeasible. 


is the computational complexity for solving the information 
subproblem. 


IV. Optimal Power Allocation with Fixed 
Eigenchannel Assignment 

In this section, we develop a waterfilling-like procedure to 
solve the optimal power allocation problem for a given eigen¬ 
channel assignment. Addressing this problem is essential for 
the development of a low complexity algorithm for problem 
©, discussed in Section [Vl it is also used to find the optimal 
power allocation of the MILP approach upon derivation of the 
eigenchannel assignment. 

Let, X and f, denote the sets of eigenchannels assigned for 
information and EH respectively, and Pf denote a fixed value 
for the total power that appears in constraint (l5bl) : then, the 
optimal power allocation problem can be defined as: 


min ^ Pi (17a) 

ie{xus} 

S>«g(l + T^j^)s«, (17b) 

iex ^ ' e J / 

Pj) > Ceh, ( 17 c ) 

iG£: je{xus} 

0 < Pi < Pmax, i G {X U £}. (17d) 


In formulation (fTTl) . the information constraint (I17bl) is 
independent of the EH constraint as no power allocated to 
energy eigenchannels appears in it, while the EH constraint 
only depends on the total power allocated to information 
eigenchannels so that the two subproblems can be solved 
almost independently. 

Note that by fixing the total power in the information con¬ 
straint of the power allocation problem, a bisection procedure 
has to be followed to obtain the optimal total power. This 
procedure is summarized in algorithm [3 where the solution 
of the information and energy subproblems is discussed in 
sections ITV-Al and |IV-B[ respectively. The computational com¬ 
plexity of algorithmic is equal to 0{\og{l/e)N\ogN), where 
log(l/e) is the estimated number of iterations of the bisection 
procedure with a stopping tolerance e, while 0{N log N) 


A. Information subproblem solution 

For given assignment and fixed total power in the informa¬ 
tion constraint, the information subproblem is expressed as: 


min Pi 
iex 


;.t. J2log(l + j 

iex ^ 

0 < Pi < Pmax 


Pi.Xi 


>Cu 


i G X. 


(18a) 

(18b) 

(18c) 


The solution of problem (IT^ is summarized in theorem 2. 


Theorem 2. Let Oi = (1 + crgPj)/A^. The solution of 
the information subproblem dll, if the problem is feasible 
log(l -|- Pmax/^i) ^ Cl), A given by. 


P* = 


0 , < Oi 

^ Oi^ Oi X: ^ — Oi ~{~ Pmax 
Pmax, Otherwise, 


where the Lagrange multiplier v is derived from: 


V = and p = 


TlieX24 + PmaxUi) 


(19) 


( 20 ) 


Sets Xi and X 2 are defined as Xi = {i \ ^ < P* < Pmax} 
and X 2 = {i: Pi = Pmax}- 


Proof: The proof is given in Appendix A. ■ 

Theorem 2 indicates that the computation of the optimal 
power allocation can be easily obtained from (O, provided 
that the value of u satisfying (l20l) is found. This can be 
achieved if the eigenchannels in sets Xi and X 2 are found. 
Towards this direction, first we need to sort the 2|X| values 
Oi and Oi -i- Pmax in ascending order; let us assume that this 
order is 9^1 < q2 <---< q2\x\- ^ ^ ^^ 1 ^ Pi nnd X 2 

can be constructed and a candidate value, z/c, can be computed 
from (I 20 I) : in case Uc G [qj,qj-^i] then it is optimal, otherwise 
the subsequent range of values needs to be examined, i.e., 
[Qj-\-i,Qj-\-2]- For example, if |X| = 2, Oi = 1, O 2 = 2, 
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Pmax = 0-5, and u = 2.3, with qi = {1,1.5, 2, 2.5}, then, 
according to (l20b , it will be true that Xi = {2} and X 2 = {1}, 
as U 0\ > Pmax and 0 <C z/ O 2 Pmax- 

A computationally efficient method to obtain the optimal 
value of u is to sequentially keep updating p and r = |Xi | for 
cheap computation of Uc, while examining new regions. When 
examining region [qj^qjj^i], two cases can occur: (a) addition 
of some eigenchannel, say i, into set Xi (when qj = Oi), 
and (b) addition of eigenchannel i into set X 2 (when qj = 

+ Pmax)^ The first case, requires addition of element i into 
set Xi so that r has to be increased by one, and p has to be 
multiplied by qj = Oi. The second case, requires deletion of 
eigenchannel i from set Xi and addition into X 2 so that r has 
to be decreased by one; updating p requires division by Oi, due 
to the deletion of eigenchannel i from set Xi , and division by 
(1 + Pmaxl^i)^ due to addition of element i into X 2 , yielding 
division by 0i(l^Pmax/0i) = Oi+Pmax = Qj - To summarize, 
starting from p = 2^^ and r = 0 the two cases require setting 
r = r + 1, p = pqj, and r = r — 1, p = p/qj, respectively; 

once r and p are updated, we compute Uc = and check to 

see whether z/c G in which case we stop, otherwise 

we examine the subsequent range of values. The complexity 
of solving the information subproblem is 0{N \og{N)), due 
to the sorting of the qi values. 

B. EH subproblem solution 

Once the information subproblem is solved, the optimal 
solution of (Ell is obtained by solving the EH subproblem 

min E Pi (21a) 

P Pi)) P Ceh, (21b) 

^ P Pi P Pmaxi i £ (21c) 

where Pi = The solution of ([21]) falls within the 

following four cases: 

Case 1: The minimum power allocated for the information 
eigenchannels, Pj suffices to satisfy the EH constraint, i.e., 
\E\Picr‘^ > Ceh- Hence, no power needs to be allocated to 
EH eigenchannels. 

Case 2: The minimum power allocated for the informa¬ 
tion eigenchannels and the total power allocated for the 
EH eigenchannels suffices to satisfy the EH constraint, i.e., 
Pmax Eies + \^\Pmax)cr'i > Ceh- Because the 

EH constraint (I17cl) is linear, power allocation is performed 
by sequentially filling EH eigenchannels, starting from the 
one with the largest eigenvalue to the one with the smallest 
eigenvalue, until the EH constraint is satisfied. 

Case 3: The total power allocated to the information 
and EH eigenchannels suffices to satisfy the constraints, 
i.e., \£\NPmax P Pmax ^ies — Ceh- In this case 
all EH eigenchannels are allocated maximum power, i.e., 
Pi 2 = Pmax, h ^ E, and extra power is arbitrarily added 
to information channels until the EH constraint is satisfied. 

Case 4: The total power allocated to the information and EH 
eigenchannels is not sufficient to satisfy either the information 
or the EH constraint, hence the problem is infeasible. 


Algorithm 4 : Group Eigenchannel Assignment 

1: Init. Set P = 00. 

2: Compute i e JV, using algorithmic 

3: if a\\{Ppf < Praax) then 

4: p* = £ = £^^I 

5: else 

6: for i e Af do 

7: Compute the minimum number of information eigen¬ 
channels with X = {i, ..., N}. 

8: for j = ip - 1,..., N do 

9: Set Xc = 

10: Set Ec = {1, - 1, j + 1,AT}. 

11: Obtain a lower bound to the total power, P^^ = 

{y — l)/crg, where y is computed from (O ac¬ 
cording to theorem 1. 

12: if (Pi^J < P) then 

13: Solve problem (fTTl) to obtain G M. 

14: if (5]]. PP^^ < P) then 

15: Store current solution: P = J]]- P^^^\ 

^ pOPA^ 


Note that the above four cases are examined sequentially 
until one of them is satisfied, yielding the optimal solution to 
dni). The complexity of solving the EH subproblem is 0{N). 

V. Low-complexity heuristic 

In general, the solution methods proposed for problem ^ 
in Sections IIII-AI and IIII-BI have exponential complexity due 
to the nonlinear and combinatorial nature of the problem, 
and hence not suitable for real-time execution. In this section 
we propose a polynomial complexity heuristic that provides 
suboptimal results, but is suitable for real-time execution. 

The heuristic is based on the observation that one con¬ 
straint usually dominates over the other, acquiring the “best” 
eigenchannels, i.e., those with the largest eigenvalues. If the 
information constraint is dominant, the first k consecutive 
eigenchannels will be assigned for information decoding and 
the rest for EH, where k is a. value to be determined. 
Nonetheless, if the EH constraint is dominant, then apart from 
the “best” eigenchannels, the “worse” eigenchannels could be 
given for EH, with power only allocated in the best ones. The 
reason is that empty eigenchannels (Pi = 0) are only useful 
for EH. Note that in both cases, information eigenchannels 
remain grouped, and for this reason we call the heuristic group 
eigenchannel assignment (GEA). 

GEA heuristic is outlined in algorithm [4j Eirst, algorithm 
[2] is examined as a way to optimally solve cases of low 
transmitted power, i.e. with Pt < Pmax, i e Af. If that 
is not the case, GEA implements the main idea of assign¬ 
ing information eigenchannels in a group, as well as two 
strategies that significantly reduce the total number of power 
allocation problems that need to be solved. The first strategy 
estimates the minimum number of eigenchannels required to 
satisfy the information constraint, given that i is the 

largest information eigenchannel. is the first number 
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that satisfies the information constraint with g‘1 = 0, i.e., 

YZi'=i + Pmax^m) > Cj. The second strategy solves 

the fixed assignment problem for Pmax = oo to obtain a lower 
bound to the total power, , according to theorem 1, so that 
the standard fixed assignment problem is solved only when the 
currently best solution P is less than P^J . Note that if any 
subproblem is infeasible the returned total power is equal to 
infinity. 

GEA involves N{N — l)/2 iterations to account all con¬ 
sidered eigenchannel combinations, while the maximum com¬ 
plexity for each assignment is 0(A^ log(A^) log(l/e)) which 
is the complexity to solve (fTTl) . Note that neither of the intro¬ 
duced strategies increase the complexity of the algorithm as 
the complexity of solving (fTTl) is larger. Hence, the computa¬ 
tional complexity of the algorithm is 0{N^ log(A^) log(l/e)). 

VI. Numerical results 

Computer simulations are carried-out in order to evaluate 
the performance of the developed algorithms. Since the scope 
of the paper is to highlight the theoretical idea of SWIFT 
in point-to-point MIMO setup, only small scale fading is 
considered for the wireless medium. Specifically, for each 
different experimental configuration, the results are averaged 
over several randomly generated instances of Nr x Nr channel 
matrix in which the entries of the channel are independent 
and identically distributed ZMCSCG random variables with 
unit variance. Thus, the information rate is measured in Bits 
Per Channel Use (BPCU) and power is measured in Power 
Units (PUs). Mathematical modeling of MISOCP and MILP 
formulations was done using the Gurobi optimization solver 

ED. 

A. The effect of harvesting and information constraints 

To evaluate the effect of harvesting and information con¬ 
straints, the optimal MISOCP solution in Q is invoked to 
solve several instances of the original problem under different 
configuration setups. In all instances, an 8 x 8 MIMO channel 
is considered (thus N = S) and the energy conversion 
efficiency factor is set to S = 0.3. In each figure, the per¬ 
formance is shown in a relative manner by using the measure 
of the Normalized Average Transmission Power, which for 
each setup is defined as the ratio between the average of the 
objective value (of this specific setup) to the maximum average 
value of the objective over all the setups which are illustrated 
in the figure. 

In Fig. |3(a)| and Fig. |3(b)[ we plot the average normalized 
transmitted power of MISOCP versus the RF EH constraint 
(Ceh) and information rate constraint (Cj), respectively, 
when cTg = 0.1 and Pmax = 2 PUs. Each point in the 
figures is the normalized average value of the total power 
consumption calculated over 1000 problem instances. As 
expected, the required transmitted power increases as Ceh 
and/or Cj take higher values because the constraints are harder 
to satisfy. One interesting observation that was extracted by 
analyzing MISOCP’s eigenchannel assignment output, is that 
a significant portion of energy constraint Ceh is covered from 


information eigenchannels by exploiting the interference term 
in (O. 

In Fig. |3(c)| we plot the impact of the peak power constraint 
Pmax on the performance of MISOCP for different RF EH 
thresholds and Cj = 10 BPCU and = 0.1. As the 
bound Pmax increases, the inherent spatial diversity of the 
system is exploited in a more efficient way by the optimal 
power allocation policy. Thus, the performance of the system 
is improved, despite the fact that potentially more power is 
available to the system. For any given Ceh and Cj value, 
there is a point where any further increase in Pmax does not 
affect the performance of the system. This is because there is 
no eigenchannel in which the upper power bound constraint is 
active, e.g for all eigenchannels the power allocated amount 
is strictly less than Pmax- Clearly, this phenomenon emerges 
for larger values of Pmax as the constraints become harder. 

In Fig. |3(d)| we show the impact of channel estimation error 
cTg on the transmitted power for different RF EH thresholds 
Ceh = {2,10,20}. As can be seen the value of signifi¬ 
cantly affects the total power consumption. For a fixed value 
of Ceh, a higher estimation error means that more power 
is required to guarantee the information rate constraint. The 
increase becomes more intense for higher values of the energy 
constraint. 

B. Comparison between MISOCP and MILP 

As it was concluded by extensive simulation results, the 
performance of MIFP is perfectly matched to MISOCP in 
terms of average transmission power. This is an interesting 
result since the two approaches aim to solve the original 
problem from two completely different viewpoints. However, 
one important issue that should be taken into consideration 
is the execution time of each solution algorithm. In Fig. 
m the execution time of MIFP and MISOCP with respect 
to the size of the problem is shown, e.g the number of 
eigenchannels, for two different setup configurations. The 
bound for the approximation error in MIFP was set equal 
to Cmax = 10“^. As cau be seen, the execution of MILP is 
more robust compared to the execution time of MISOCP which 
becomes non-practical as N increases. Having algorithms that 
provide robust execution times for large N, such as the MILP 
algorithm, is particularly important in obtaining the optimal 
solution in massive MIMO configurations where hundreds of 
antenna elements are present. 

C. Comparison between MISOCP and GEA 

In Fig. |5(a)[ the performance of the Group Eigenchannel 
Assignment (GEA) algorithm and MISOCP is compared by 
averaging 10000 problem instances. An 8x8 MIMO channel is 
considered with Pmax = 2 PUs, = 0.1 and several different 
values Ceh and Cj. As can be seen, GEA has excellent 
performance, achieving results within 5% from the optimal 
MISOCP solutions for all the examined configuration setups. 
This is a quite significant result if we take into consideration 
that GEA is a low-complexity algorithm, while MISOCP 
has exponential complexity. The small gap in performance 
is because GEA cannot output an eigenchannel assignment 
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where information eigenchannels are not sequential, which 
appear to be the optimal assignment in certain cases. In 
Fig. |5(b)[ the feasibility of the two algorithms is compared 
against the Ceh parameter for three different values of Cj. 
Again, a small gap exists in GEA’s feasibility which becomes 
more notable when the RF EH and/or information constraints 
become more stringent. This observation is inline with the 
results in Fig. |5(a)| where the performance of GEA deteriorates 
as Ceh and Cj take higher values. 

VII. Conclusions 

In this paper, we have theoretically investigated SWIFT in 
the spatial domain for a MIMO channel with RF EH capabil¬ 
ities. By using SVD decomposition for the wireless channel 
under an imperfect channel knowledge, the proposed technique 
uses the eigenchannels to convey either information or energy, 
with the goal of minimizing the overall transmitted power 
subject to information and energy constraints. Although the 
examined problem is non-linear and combinatorial, MISOCP 


and MILP formulations have been developed that provide 
optimal and very close to optimal solutions, respectively. 
Using Lagrange theory, a waterfilling-like procedure for the 
optimal power allocation is derived when the eigenchannel 
assignment is known. Finally, a polynomial complexity algo¬ 
rithm is presented which produces a near-optimal solution for a 
wide range of parameter configurations. In this paper, we have 
laid the theoretical foundation of the spatial domain SWIFT 
which can be a promising technology especially for millimeter 
wave communication systems. In future work, practical imple¬ 
mentations of such systems need to be investigated. 

Appendix A: Proof of Theorem 2 

In order to prove the theorem, we consider the Karush- 
Kuhn-Tucker (KKT) optimality conditions, as (fTSl) is a convex 
optimization problem ll40ll . Let /if and /if denote the La¬ 
grange multipliers corresponding to the constraints Pi > 0 and 
Pi < Pmax^ i G X, respectively and u the Lagrange multiplier 
corresponding to the information constraint (IlSbl) . Then, the 
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Fig. 4. Execution time versus number of eigenchannels for different RF EH 
and data constraints. 


KKT conditions can be written as: 



(22a) 

iex 

0<Pi< Pmax, i e I, 

/ \ 

(22b) 

i/(Vlog(l + P,/0i)-C/1 =0, 

(22c) 

\iex / 

/if Pi = 0,i G X, 

(22d) 

{Pi — Pmax) = 0, i G X, 

(22e) 

z/ > 0, /if > 0, /if > 0, i G X, 

(22f) 

1 — z// (Oi + Pi) — /if + /if = 0, i G 

(22g) 


restrictions hold: 

Restriction 1: Because Cj > 0, the minimization of 
implies that constraint (I22al) is always satisfied with equality 
at the optimal solution; otherwise, some Pi could be further 
reduced to achieve equality for the information constraint 
and also reduce the objective function value. Hence, (I22cl) is 
always satisfied and > 0. 

Restriction 2: Eqs. (I22dl) and (I22eb imply that when = 0 
it is true that = 0, when 0 < P^ < Pmax it is true that 

/if = 0 and /if = 0, and when Pi = Pmax it is true that 

nY = 0 . 

Restriction 3: According to the values of /if and /if we can 
identify three cases for the value of u obtained from ( |22g| ): 

3.1 Case /if = 0 and /if > 0 (Pi = 0): /if = 1 — vjOi > 

^ u < Oi. 

3.2 Case = /if = 0 (0 < P^ < PmaxY + P^ > 0. 

3.3 Case /if = 0 and /if > 0 (P^ = PmaxY > Oi P 
Pmax P 0- 

Restriction 3, guarantees that (O is satisfied for all infor¬ 
mation eigenchannels and that zi > 0 if P^ > 0, for some 
i G X, which is always true since (7/ > 0. 

Notice that the contribution of eigenchannels with P^ = 0 
or Pi = Pmax is not depended on u (it is equal to 0 and 
\og(lPPmax/Oi), respectively), while for 0 < P^ < Pmax, the 
contribution is equal to \og{u/Oi) (obtained from substitution 
of Pi = u — Oi into the associated log term of the information 
constraint). Substituting all contributing channels into (I22al) 
yields: 

y] l0g{iy/di) + l0g(l + P^a:c/0i) = Cl 

where Xi and X 2 are defined in Theorem 2. Rewriting the 
logarithmic terms in product form yields: 


where (I22al) - (l22bb are the primal feasibility constraints, (I22cb - 
(I22eb are the complementary slackness conditions, (I22fb are 
the dual feasibility constraints, while ( |22g| ) is obtained by 
setting the derivative of the Lagrange function with respect 
to Pi to zero. Based on the KKT conditions ([22b a number of 


f 1-r I 

log { =-- n (1 + Pma./ei) ^ = Cl 

which after simple algebra yields (l20b . This completes the 
proof. 
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